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Abstract 

We calculate branching ratios for pure penguin decay modes, B — ► ^if decays using perturbative 
QCD approach. Our results of branching ratios are consistent with the experimental data and larger 
than those obtained from the naive factorization assumption and the QCD-improved factorization 
approach. This is due to a dynamical penguin enhancement in perturbative QCD approach. 
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1 Introduction 

Recently the branching ratios of B — ► <pK decays have been measured by the BaBar Q|, BELLE g 
and CLEO || collaborations. There is an interesting problem related a penguin contribution to decay 
amplitudes B. A naive estimate of the loop diagram leads to P/T ~ a s /{12i:)\og{m1/m 2 c ) ~ 0(0.01) 
where P is a penguin amplitude and T is a tree amplitude. But experimental data for Br(S — » ifor) and 
Br(J3 — ► 7T7r) leads to P/T ~ O(0. 1). Therefore, there must be a dynamical enhancement of the penguin 
amplitude. This problem is studied by Keum, Li and Sanda using perturbative QCD (PQCD) approach 
0.5^ <f)K modes are important understanding penguin dynamics, because these modes are dominated 
by penguins. Here we report our study of B — > <f>K decays using PQCD. 

PQCD method for inclusive decays was developed by many authors over many years, and this for- 
malism has been successful. Recently, PQCD has been applied to exclusive B meson decays, B — > Kir 
|j||, 7T7T pj, 7rp, -kuj J7p, KK H| and Krf- > 0. PQCD approach is based on the three scale factorization 
theorem [[L0|, pjj|. For example, B —* K transition form factor can be written as 



? BK 



[dx)[db}Ci(t)$ K (x2,b 2 )H(t)$B(xi,h)exp 



- V / — t0(« s (m)) 



(i) 



where X\ and x 2 are momentum fractions of partons, b\ and b 2 are conjugate variables of parton transverse 
momenta hyr and k 2 T, and 70 is the anomalous dimension of mesons. The hard part H(t) can be 
calculated perturbatively. d (t) is the Wilson coefficient corresponding to the four-quark operator causing 
B — > K transition. The scale t is given explicitly in terms of x\, x 2 , b\, b 2 and Mb, and it is of O(vAMb). 
Here A = Mb — rnb, where Mb and nib are B meson mass and b quark mass, respectively. It is important 
to note that in PQCD, the scale of the Wilson coefficient t can reach below Mb/2. In the factorization 
assumption p2[ , this scale is fixed at Mb/2. The Wilson coefficient for a penguin operator increases as 
the scale evolves down. This explains the enhancement of the penguin amplitude in PQCD compared to 
the amplitude obtained by the factorization assumption. 
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In this method, we can calculate not only factorizable amplitudes but also nonfactorizable and annihi- 
lation amplitudes. In case of B — ► <pK decays, the factorizable amplitudes which can be written in terms 
of form factors F BK and F^ K are shown in Fig. fi](a)-(d). The nonfactorizable amplitudes are shown in 
Fig. [|(a)- (d). Ellipses denote meson wave functions in these figures. For illustration purposes, we show 
the hard part of the nonfactorizable diagram as the dashed box in Fig. 0(a). The parameters in meson 
wave functions are calculated from the light-cone QCD sum rules, and the theoretical uncertainty of the 
parameters is about 30%. The hard part depends on the particular processes, but it is calculable. The 
wave functions contain non-perturbative dynamics and are not calculable, but once it is known, it can 
be used for other decay processes. 

In this paper, we calculate branching ratios for B — ► <fiK modes using PQCD approach. The detail 
is discussed in Ref. [[l3|. We predict the branching ratios for B — > <\>K decays, and our predictions agree 
with the current experimental data and are larger than the values obtained from the naive factorization 
assumption (FA) and the QCD-improved factorization (QCDF) Q, |lo| . 
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Figure 1: Feynman diagrams contributing to factorizable amplitudes for B — > cf>K 
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Figure 2: Feynman diagrams contributing to nonfactorizable amplitudes for B — > <f>K 



2 B — > cj)K Amplitudes 



We consider B meson to be at rest. In the light-cone coordinate, the B meson momentum Pi, the K 
meson momentum P2 and the <fi meson momentum P3 are taken to be 
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Pi = -£(1,1, Or), P 2 = -^-(l-H,0,0 T ), P 3 = -^(ri,l,0 T ), 



V2 



V2 



V2 



(2) 



where r^ = M^/Mb, and the K meson mass is neglected. The momentum of the spectator quark in the 
B meson is written as ki. Since the hard part is independent of k^ , the S(k^~) function appears after 



integrating over its conjugate spacial variable. Therefore, k\ has only the minus component fcj" and small 
transverse components ki-p. k^ is given as fcj" = x\P^ , where x± is a momentum fraction. The quarks in 
the K meson have plus components x 2 P 2 and (1 — ^2)^2 > and the small transverse components k2T and 
— \i.2T, respectively. The quarks in the <fi meson have the minus components X3P3 and (1 — Xs)P^ , and the 
small transverse components k 3 T and — k 3 T, respectively. The <p meson longitudinal polarization vector e^ 
and two transverse polarization vector e^t are given by e^ = (l/-\/2?>)(— r\, 1, Ot) and e^t = (0, 0, It). 
The B meson wave function for incoming state and the K and <fi meson wave functions for outgoing 
state with up to twist-3 terms are written as 



^B,af3,ij 



J d^kiTe-^ 1 ^- 1 ^ 21 ^ [(A + M s ) 75 <Mxi,k 1T )] a/3 » (3) 

V zi\ c Jo 
^L = "S= f dx^ p ^ [M+ UM^)+ h /P 3 <^(* 3 ) + M^(* 3 )] , (5) 

V^IVc Jo ^ 

where i and j is color indices, and a and /3 are Dirac indices, vhqk is related to the chiral symmetry 
breaking scale, ttiqk = Mj|-/(md+ms). w and n are defined as v^ = P 2 jP 2 and n M = z/j 1 1 z 2 = (0> 1> Ot)- 
We neglect the terms which are proportional to the transverse polarization vector eT, because these terms 
drop out from our calculation kinematically. The explicit form of these wave functions will be shown in 
Sec. I 

Widths of B — ► <j)K decays can be expressed as 

^- ° F |-4| 2 . (6) 
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The decay amplitudes, A, and A, corresponding to B° — > <pK°, and B° — > </>-K"°, respectively, are written 

as 

A = uv ts v; b F e p + v ts v; b M P + f B V ts V t * b F p + v ts v; b M p , (7) 

A = UV t * s V tb F e p + V t * s V tb M P + f B V t * s V tb F p + V t * 3 V tb M p , (8) 

where F e is the amplitude for factorizable diagrams which are considered in FA. F a and M. are the 
annihilation factorizable and the nonfactorizable diagrams which are neglected in FA. The indices, e, 
and a, denote the tree topology, and annihilation topology, respectively. The index P denotes the 
contribution from diagrams with a penguin operator. The decay amplitudes, A + , and A~ , corresponding 
to B + — > 4>K + , and B~ — > 4>K~ , respectively, are written as 

a + = uv ts v; b F P + v ts v; b M P + f B v ts v t iF p + v ts v; b M p - f B v us v; b F? - v us v: b M T a , (9) 
a- = UV t *v tb F p + v; s v tb M P + f B v; s v tb F p + v; s v tb M p - f B v; s v ub F^ - v: s v ub M T a , (io) 

where the index T denotes tree contributions. Since the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements for the tree amplitudes are much smaller than for the penguin amplitudes, the tree contributions 
are very small. 

The factorizable diagrams are given as Fig. nl The factorizable penguin amplitude, F p , which comes 
from Fig. |l|(a) and Fig. |](b) is written as 

F p = 8nCpMg / dx\dx 2 \ b\db\b 2 db 2 4> B {x\,b\) 

Ja Jo 

x ( [(1 + 12)^(12) + r K (l - 2x 2 ) {(feixi) + <&fa))] E e (tW)N t {x 2 {l - x 2 )} c h e (x 1 ,x 2 ,b 1 ,b 2 ) 



+2r K ^(x 2 )E e (t^)N t {x 1 (l - x 1 )} c h e {x 2 ,Xub 2 , h) I , 



(11) 

where N t {x(l — x)} c is the factor for the threshold resummation jl(|. We use N t — 1.775 and c = 0.3 
|p7[ . The evolution factors are defined by E e (t) = a s (t)a e (t) exp[— 5e(i) — S'if(t)] where exp[— S'i(i)] 
is the factor for the kx resummation flq], [ J19J . The explicit forms of the factor Si(t) are given, for 
example, in Ref. ||. The hard scales t e and t e , which are the scales in hard process, are given by 
V e ' — max(y / X2Ms, l/i»i, I/&2) and 4 J = max( v / xiAfs, I/61, 1/62)- The Wilson coefficient is given by 






(12) 



The hard function, which is the Fourier transformation of the virtual quark propagator and the hard 
gluon propagator, is given by 



h e (x u x 2 ,b ll b 2 ) = K (^xTx^M B b 1 )[9(b 1 ~b 2 )K (^M B b 1 )I (^M B b 2 ) + {b 1 ^b 2 )] .(13) 

The factorizable annihilation diagrams shown in Fig. [j](c) and Fig. 0(d), and the nonfactorizable 
diagrams shown in Fig. @(a)-(d) can be also calculated in the same way as F^ fll3||. 



3 Numerical Results 



We use the model of the B meson wave function written as 
4>B(x,b) = NbX 2 (1 — x) 2 exp 



1 (xMbV uj B b 2 
2\ lub J 2 



where lob = 0.40 GcV p0| . Nb is determined by normalization condition given by 
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The K meson wave functions are given as 
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<(> T K {x) = 
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1 + 6 ( 5?73 - - 773^3 - 20 Pk - g P/f ft 2 ) (1 10* + I0x z ) 



(14) 
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where px = {m-d + ins) I Mr pi| ) p2| - The parameters of these wave functions are given as ai 
0.17, a 2 = 0.20, 773 = 0.015 and cj 3 = —3.0 where the renormalization scale is 1 GeV. 
The <j) meson wave functions are given as 
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<l><t>(x) 
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2J2K. 



:6a;(l — x) , 
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(20) 
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where £j = 0.024 and 5+ = 0.46 |23|] . We have found that the final results are insensitive to the values 
chosen for £j and 5+. 

We use the Wolfenstein parameters for the CKM matrix elements A = 0.819, A = 0.2196, i?& = 
\J p 1 + rf = 0.38 |24|], and choose the angle 4>z = n/2 [g|. We have found that the final results are quite 
insensitive to the values of ^3. For the values of meson masses, we use Mb = 5.28 GeV, Mk — 0.49 GeV 
and M^ — 1.02 GeV. In addition, for the values of meson decay constants, we use fs = 190 MeV, fx — 
160 MeV, fy = 237 MeV and fj = 215 MeV. The B meson life times are given as t b o = 1.55 x 10~ 12 sec 
and t b ± = 1.65 x 10~ 12 sec. And we use A^ D = 0.250 GeV and m 0K = 1.70 GeV §. 

We show the numerical results of each amplitude for B° — > cj}K° and B ± —* <\>K^ decays in Tab. |[ The 
factorizable penguin amplitude F^ gives a dominant contribution to B —> <j)K decays. The factorizable 
annihilation penguin amplitude F^ generates a large strong phase. In B — * cf>K modes, there are 
contributions from fsFj and Mj . These tree amplitudes contribute only a few percent to the whole 
amplitude, since the CKM matrix elements related to the tree amplitudes are very small. In order to 
isolate the trivial uncertainty from fs, fa and f$, we express our prediction for B — > cf>K as 



Bi{B° -> 4>K°) = 
Br(B ± -» c^K^) = 



fBfxfd) 



190 MeV 160 MeV 237 MeV 



x (9.43 x 10~ 6 ) , (22) 



x (10.1 x 10~ 6 ) . (23) 



190 MeV 160 MeV 237 MeV 

We found that our result is insensitive to fT /fy- For example, 10% variation of fT/fy leads to less than 
1% variation in our final result. The current experimental values are summarized in Tab. g. The values 
which are predicted in PQCD are consistent with the current experimental data. However, our branching 
ratios have the theoretical error from the 0(a 2 s ) corrections, the higher twist corrections, and the error of 
input parameters. Large uncertainties come from the meson decay constants, the shape parameter u>b, 
and itiqj{. These parameters arc fixed from the other modes {B — » Kn, Dn, etc.). We try to vary u>b from 
0.36 GeV to 0.44 GeV, then we obtain Br(B ± -> (pK^ = (7.54 - 13.9) x 10~ 6 . Next, we set uj b = 0.40 
and try to vary m 0K from 1.40 GeV to 1.80 GeV, then we obtain Br(B ± -> (pK^) = (6.65 - 11.4) x 10~ 6 . 
Now, we consider the ratio of the branching ratio for the B° — > <pK° decay to the one for the 
B + — » 4>K + decay. The theoretical uncertainties from various parameters are small, since the parameters 
in the numerator cancel out those in the denominator. The difference between the two branching ratios 
come in principle from B meson life times, tree and electroweak penguin contributions in annihilation 
amplitudes. We found that the tree and electroweak penguin amplitudes in the annihilation diagrams are 
negligible. Tree amplitudes are suppressed by two factors. First, they are annihilation processes which 
are suppressed by helicity. Second, they are multiplied by small CKM matrix elements. We predict that 
this ratio is 

Br(£° -> <f>K°) 



Bi(B+ -c 4>K+) 



= 0.95 , (24) 



where the theoretical uncertainties from m^K and lub are less than 1%. The ratio is essentially given by 
the life time difference. The experimental value of this ratio from BELLE g] is Br(i3° — ► 4>K°)/Br(B + — > 
<f)K+) = 0.82±°;i|±0.10. 

In FA, the branching ratio is very sensitive to the effective number of colors N*ff. If we set N** ■* =3, 
then the branching ratio is about 4.5 x 10~ 6 where the scale of the Wilson coefficient is taken to Mg/2 
and F BK is 0.38 from the BSW model. In QCDF, branching ratios for B — * <pK decays are predicted 
as Bt(B° -c (j)K Q ) = (4.0+ 2 ;^) x 10~ 6 and Bi(B~ -> ^>K~) = (4.3l?;2) x 10_6 including the annihilation 
diagram p5[. Our predicted values are larger than these results. This is due to the enhancement of the 



Wilson coefficient for the penguin amplitude as explained in Sec. |j. In PQCD approach, the scale of the 
Wilson coefficients, which is equal to the hard scale £, can reach lower values than Mb/2. 

4 Summary 

In this paper, we calculate B° — > (f>K° and B^ — > ipK^ decays in PQCD approach. Our predicted 
branching ratios agree with the current experimental data and are larger than the values obtained by FA 
and QCDF. Because the Wilson coefficients for penguin operators are enhanced dynamically in PQCD. 

Note added: 

After this work has been completed, we become aware of a similar calculation by Chen, Keum and Li 
p5| . Our results are in agreement. 
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B" -v <j)K° B ± -> QK* 



f^Ff -1.03 x 10- 1 -1.03 x 10- 1 

f B F? 6.45 x 10~ 3 + i 4.28 x 10~ 2 6.17 x 10~ 3 + i 4.20 x 10~ 2 

M e p 5.24 x 10~ 3 - i 3.61 x 10~ 3 5.24 x 10~ 3 - i 3.61 x 10~ 3 

Mf -8.03 x 10~ 4 - i 1.73 x 10~ 3 -6.56 x 10~ 4 - i 7.22 x 10" 

j B Fl -1.11 x 10" 1 - i 3.75 x 10- 

Ml 1.60 x to" 2 + i 2.77 x 10" 2 



Table 1: Contribution to B° — > <j)K° and S ± — > §K ± decays from each amplitude 





Br(B -> (f>K°) 


Br(B ± -> QK^) 


BaBar 


(8.1±l;g =b 0.8) x 10" 6 


(7.7t\:t ± °- 8 ) x 10 " 6 


BELLE 


(8.7±l* ± 1.5) x 10- 6 


(10.6l 2 ;g±2.2) x 10~ 6 


CLEO 


< 12.3 x 10~ 6 


(5.5l 2 ;g±0.6) x 10" 6 



Table 2: The experimental data of B -* <f>K branching ratios from BaBar§, BELLE@ and CLEO(| 



